Quark mean field model for nuclear matter and finite nuclei 
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Abstract 

We study nuclear matter and finite nuclei in terms of the quark mean field (QMF) 
model, in which we describe the nucleon using the constituent quark model. The meson 
mean fields, in particular the a meson, created by other nucleons act on quarks inside a 
nucleon and change the nucleon properties in nuclear medium. The QMF model predicts 



X 

an increasing size of the nucleon as well as a reduction of the nucleon mass in the nuclear 
environment. The present model is applied to study the properties of finite nuclei after 
fixing all the parameters by the nuclear matter properties, and it is found to give satis- 
factory results on the nuclear properties. 
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1 Introduction 



We have been describing the nucleus with the assumption that the nucleon properties 
are unchanged from those of the free space nucleon. However, after the EMC effect 
was reported, we encounter often the discussions on the change of hadron properties in 
nuclei (lj . The EMC effect motivated many theoretical works on the study of hadrons in 
terms of quarks and gluons, although the pion interpretation is not excluded ||]. Among 
them, the QCD sum rule suggested a large change of the masses of the vector mesons and 
nucleons in nuclei ||. It is very important to pursue these possibilities both theoretically 
and experimentally from various view points. 

On the other hand, recent theoretical studies show that the properties of nuclear mat- 
ter can be described nicely in terms of the relativistic Brueckner-Hartree-Fock (RBHF) 
approach With the nucleon- nucleon interaction fixed from the nucleon- nucleon scat- 
tering and the deuteron properties, we can reproduce very closely the nuclear matter 
saturation properties. The reason of the success is the large reduction of the nucleon 
effective mass in nuclear matter, which provides a density dependent repulsive contribu- 
tion for the energy density. This reduction of the nucleon mass is also the source of the 
large spin-orbit splitting in finite nuclei, which is the key phenomenon of the nuclear shell 
model. Any model which changes the hadron properties in nuclear matter should respect 
these observations. 

Guichon proposed an interesting model on the change of the nucleon properties in 
nuclear matter [fJ. The model construction mimics the relativistic mean field theory, 
where the scalar and the vector meson fields couple not with nucleons but directly with 
the quarks in nucleons of the nuclei ||. Hence, the nucleon properties change according 
to the strengths of the mean fields acting on the quarks. The nucleon is modeled in 
terms of the MIT bag model |7]]. The scalar meson provides a strong attraction and 
as a consequence provides a negative mass to the quarks. The MIT bag model then 
reduces the nucleon mass in nuclear matter. This model was extended by Thomas and 
his collaborators under the name of the quark-meson coupling (QMC) model and applied 
to many observables || . 

The model of Guichon relies strongly on the choice of the nucleon model, for which 
the QMC model uses the MIT bag approach Q. The MIT bag model assumes that the 
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inside of the bag is the perturbative vacuum and the quark mass is the bare mass, which 
is nearly zero for the up and down quarks. The perturbative vacuum has a larger energy 
than the non-perturbative vacuum and therefore the bag constant takes care of the energy 
difference between the two vacua. The non-perturbative objects as mesons are not allowed 
to exist in the bag interior and stay outside the bag. The chiral symmetry, the continuity 
of the axial current, requires the coupling of quarks with pions at the bag surface |J. 

There are several conceptual problems to take the MIT bag model in the Guichon 
model; 

1. The non-perturbative objects, a and lo mesons, have to be present in the perturbative 
vacuum. 

2. The quark mass has to change from its bare mass due to the coupling to the a meson. 

3. For the up and down quarks, the resulting mass of the quarks is negative. 

There may be some arguments to overcome the above mentioned points and justify the 
QMC model §. 

We would rather like to take another model for the nucleon, which is the constituent 
quark model . In this model, the quarks get constituent quark masses due to sponta- 
neous chiral symmetry breaking. It is then natural to have nearly zero mass pions as the 
Nambu-Goldstone bosons. Their coupling to the constituent quarks is provided by the 
chiral symmetry. This consideration makes a simple interpretation of the direct coupling 
of not only pions but also other mesons as a and u mesons, since there is no boundary to 
separate the perturbative and the non-perturbative regions in this picture. The a mean 
field could be considered as the amount related with the change of the chiral condensate 
in nuclear medium. Hence, it is natural to get the reduction of the quark mass in the 
nucleon inside of nuclei from the quark mass of the nucleon in the free space. 

The constituent quark model is used extensively also for nucleon-nucleon interaction 
and later extended to baryon-baryon interactions in the SU(3) sector with great suc- 



cess [11, 12, 13, 14, 15, 16|. In this picture, it is natural that the mesons couple with 
quarks, since the nucleon is the composite of quarks. 

Hence, it is very interesting to construct the Guichon model, where the nucleon is 
described in terms of constituent quarks, which couple with mesons and gluons. This 
model (we refer it as the quark mean field (QMF) model as named in Ref. [fT^] .) has a direct 
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connection to the one-boson exchange model (Bonn potential). We expect, therefore, the 
quantitative results similar to those of the RBHF theory Q . 

The paper is arranged as follows. In Section 2 we describe the concept of the QMF 
model and the procedure to perform calculations. In Section 3 we show the results for 
nuclear matter and in Section 4 the results for finite nuclei are presented. Section 5 is 
devoted to the summary of this paper. 



2 Quark mean field model 



We are still far away from describing nucleons and nuclei in terms of quarks and gluons 
using QCD. This difficulty arises from the fact that the theory becomes highly non- 
perturbative at low energy {p < lGeV). In addition, the QCD vacuum is realized in 
a non-trivial way, where the chiral symmetry is broken and the quarks and gluons are 
confined. The lattice QCD (LQCD) may handle hadrons. It is, however, too much to 
ask the LQCD to describe a system of many nucleons, the nuclei. Hence, we resort to an 
effective theory of QCD at low energy, which is based on QCD. The dual Ginzburg-Landau 
(DGL) theory may be a good candidate . 

In this paper, we take a more phenomenological view point. We shall begin with a 
possible Lagrangian of the quark many body system. For this Lagrangian, we take into 
account the consequence of the non-perturbative gluon dynamics of spontaneous chiral 



symmetry breaking and quark confinement [19]. The effective Lagrangian of this level 
may be written as 

C = q(i^ -m q -xc- 91^ - gffaTi) q + C( Xc , A^ . (1) 

Here, q denotes the quark fields with constituent quark mass m q , which is of the order 
of 300MeV to be consistent with the quark condensate and the pion decay constant ]p0[. 



The confinement is expressed in terms of Xc which is given by the gluon dynamics [18|. 
The path from QCD to this expression is not yet known. One method might be the use of 
the dual Ginzburg-Landau (DGL) theory, which contains QCD monopoles in the Abelian 



space in the Abelian gauge and their condensation to induce the dual Meissner effect [18|. 

are the gluon fields with the running coupling constant g at the model scale. From the 
lattice QCD data, we find the effect of the gluon exchange interaction below r < 0.2 fm and 



hence the running coupling constant should be taken at the momentum scale, p ~ lGeV. 
We include meson fields (pi, which couple with quark fields with the Dirac matrices Tj and 
the coupling strength g\. The rest C(xc, A^, (pi) denotes the confinement field, gluon and 
meson dynamics, which is not specified here explicitly. The explicit form is quite involved. 

To proceed to many body system, we take the mean field approximation for the meson 
fields. We restrict to a, u, and p mesons, which are commonly used in the relativistic mean 
field model || 21, E2|. We note that in the mean field approximation, the pion field does 



not survive due to the spin average. Since quarks are confined completely in hadrons, we 
work out the many body problem in two steps. First, we construct the nucleon under the 
influence of the meson mean fields. Then in the second step, we solve the entire nuclear 
system with the change of the nucleon properties due to the presence of the mean fields, 
which is obtained in the first step. 

The first step is to generate the nucleon system under the influence of the meson mean 
fields. In the constituent quark model, the quarks in a nucleon satisfies the following 
Dirac equation: 

-m q -Xc- 9 q Mr) - g q Mr) 7 ° - g q p p(r)T 37 °] q(r) = 0, (2) 

where t% is the isospin matrix in our nuclear physics convention. Assuming the meson 
mean fields are constant within the small nucleon volume, we can then write the Dirac 
equation as 

■ia ■ V + f3m* q + Pxc q(r) = e*q(r), (3) 

where m* = m q + g^a and e* = e — g^oj — g p pT3, with a, oj, and p being the mean fields 
at the middle of the nucleon. e is the energy of the quark under the influence of the a, 
to, and p mean fields. The quark mass is modified to m* due to the presence of the a 
mean field. Here, g q , g^,, and g q p are the coupling constants of the cr, uj, and p mesons 
with quarks, respectively. We take into account the spin correlations, E sp i n , due to gluons 
and pions so that the mass difference between A and nucleon arises. Hence, the nucleon 
energy is expressed as E* = 3e* + E sp i n , where the vector contribution is removed here. 
There exists the spurious center of mass motion, which is removed in the standard method 



by M* = \J E* 2 — < pl m >, where < p^ m >= Ya=i < Pi >> since the three constituent 
quarks are moving in the confining potential independently. 



We now move to the second step, in which the nuclear many body system will be 
solved with the change of the nucleon properties obtained in the first step. We assume 
the following QMF Lagrangian, 

£-qmf = $ - M* - g^uj - g p pT 3 j° ip + £m(o", p). (4) 

Here, ip denotes the nucleon fields. The change of the nucleon properties is worked out in 
the first step and the outcome is exclusively expressed in M*, which is a function of the 
quark mass correction, 5m q , related to the a mean field as 5m q = —g^cr. The to and p 



mean fields do not cause any change of the nucleon properties [23|, and they appear merely 
as the energy shift. These contributions are carried over as the nucleon-meson coupling 
terms with the replacement of the quark- meson couplings as g u = 3g^ and g p = 9% We 
will apply the QMF model to study the properties of nuclear matter and finite nuclei in 
the following two Sections. 



3 properties of nuclear matter 

We calculate first the change of the nucleon properties as a function of the quark mass 
correction, 5m q , which is defined as 5m q = m q — m* = —g%o. Here, the constituent 
quark mass is taken to be one third of the nucleon mass; m q = M n /3 = 313MeV. 
We take into account confinement in terms of the harmonic oscillator potential together 
with two Lorentz structures; (1) scalar potential Xc = \kr 2 (2) scalar-vector potential 
Xc = \kr 2 {l + 7°)/2. As pointed out in Ref. [p4| , the quark can not be confined when 
the vector potential is larger than the scalar one. Here, we just take two extreme types, 
since the Lorentz structure of the confinement is not established. As for the strength of 
the confining potential, we take k = 700 and lOOOMeV/ fm 2 , in order to see the results 
depending on this factor. The spin correlation, E sp i n , is fixed by the free nucleon mass 
as M n = J (3e + E sp i n ) 2 — < p 2 m > = 939MeV . We assume further that the confining 
interaction and the spin correlations do not change in the nuclear medium. 

We show in Fig.l the results of the effective nucleon mass M* as a function of 5m q . The 
results of the scalar potential with two oscillator parameters, k = 700 and WOOMeV/ fm 2 , 
are shown by the two solid curves. The effective nucleon mass decreases with bm q and 
the curvature is clearly negative, which influences the nuclear matter properties, as will 
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be discussed soon. The results of the scalar-vector potential with the same two oscillator 
parameters are shown by the dashed curves. In this case the dependence of M* on 5m q 



is almost linear as has been reported by Toki et al. [17]. We note that the results of the 
nuclear matter are very sensitive to the behavior of M*. 

To perform the nuclear matter calculation, we use the relativistic mean field (RMF) 
approximation containing nucleons, neutral scalar (c) and vector (a;) and isovector (p) 
mesons. We need to specify the meson Lagrangian Cm, since there are various versions of 
the meson Lagrangian in the RMF theory. The RMF theory with the TM1 parameter set 
includes the nonlinear terms both for a and u mesons, which can reproduce the feature of 



the RBHF theory and satisfactory properties of finite nuclei [25]. In the present model, 
we would like to take the meson Lagrangian as 

Cm = \{d,af -\mla 2 -\g z a' (5) 

1 ■■ « \ 2 i 1 2, ,2 i 1 , ,4 



Comparing with the Lagrangian in the RMF(TMl) model [25], we have deleted the nonlin- 
ear term ^g2<J 3 , so that less free parameters are contained in the present model. Actually, 
the nonlinear term \g2&^ plays similar role as the nonlinear term jg^cr 4 . Note that we 
have dropped the x p\ term, since this term vanishes at the mean field level. Now, 
the nucleons and mesons obey the Euler-Lagrange equations derived from the QMF La- 
grangian, which can be written, for the nuclear matter, as 

i^d" — M* — 9uj lo 7 - g pP r 3l °] $ = 0, (6) 
mla + g 3 a 3 = -^(W), (7) 
mluj + c 3 uj 3 =&,$7V}, (8) 

m2 P p = gpitpni ^}- (9) 

Here, the bracket ( } means the expectation value of the operator between the nuclear 
ground state. For the to and p parts we have replaced the quark-meson couplings g\ 
by the nucleon- meson couplings g% as g u = 3g^ and g p = g q p . The effective mass M* 



and its derivative with respect to the a mean field, dM*/da, are not trivial functions of 
the a mean field, because M* depends on 5m q = —g%(J non-trivially as shown in Fig.l. 
Comparing with the RMF theory, the dM*/da in the QMF model is equal to the g a in 
the RMF model. We show in Fig. 2 the quantities g a (a)/g a (0) as a function of 5m q . In 
the case of the scalar- vector potential, the dependence of g a (cr) on 5m q is very small. On 
the other hand, in the case of the scalar potential, g a ((j) increases rapidly with 5m g . 

In the present model, there are five free parameters, g%, g%, g%, C3, and g p , which 
need to be determined. We would like to follow the method in Ref. [26] to determine 
these parameters. We determine these five parameters by reproducing five equilibrium 



properties of nuclear matter [25, 26]. The equilibrium properties used here are listed 
in Table 1. As for the other parameters, we take = 783MeV and m p = 770MeV. 
We note that the variation of m a at fixed g% /m a and 53 / m 4 has no effect on the nuclear 
matter properties, but the a meson mass determines the range of the attractive interaction 
and affects the nuclear surface slope and its thickness for finite nuclei and hence the finite 
nuclear properties. The mass of the a meson is chosen to reproduce the charge radius 
of 40 Ca to be around 3.45/m. The parameter sets for the four cases used in the present 
model are given in Table 2. 

Using the parameters given in table 2, the nuclear matter properties can be obtained 
by solving the above Euler-Lagrange equations self-consistently. We plot in Fig. 3 the 
energy per nucleon, E/A, as functions of the nuclear matter density p. In Fig. 4, we plot 
the scalar potential, Us, and the vector potentials, Uy, as functions of density p. The 
results in the RMF(TMl) model are also shown for comparison. Definitely the results at 
P = Po = 0.145/m -3 do not change due to the construction of the parameter sets, but 
the results at different densities depend somewhat on the parameter set. 

Since the QMF model uses the constituent quark model to describe the nucleon in 
medium, the nucleon properties change according to the strengths of the mean fields. We 
plot the ratio of the nucleon rms radius in medium to that in free space as a function 
of density in Fig. 5. All of the four curves show the increase of the nucleon radius in 
medium. The nucleon radius increases by about 5% ~ 9% at the normal matter density 
for different confinement parameters and the type of the potential. It would be very 
interesting to calculate the EMC effect with the change of the nucleon properties, in 
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particular, with the increase of the nucleoli radius. If we make a comparison with the 
QMC model of Thomas and his collaborators II, a significantly swollen nucleon radius 
is predicted in the present QMF model without adding extra parameters, while almost 
no swelling is found in the orginal QMC model when the bag constant does not depend 
on density ||]. We note here that there are a few modified versions of the QMC model, 
which could reproduce the swollen nucleon radius by changing the bag constant p(| . 



4 Properties of finite nuclei 



To study the properties of finite nuclei, the electromagnetic field must also be included, 
which does not appear in infinite nuclear matter. The construction of the nucleon inside 
nuclei will be rather complicated if the variation of the meson mean fields over the nucleon 
volume is considered. We have to take some suitably averaged form for the meson mean 
fields in order to make the numerical solution feasible. Here, we use the local density 
approximation, which replace the meson mean fields at the quark level by their value at 
the center of the nucleon and neglect the spatial variation of the mean fields over the 
small nucleon volume. If we restrict our consideration to spherically symmetric nuclei, 
the Euler-Lagrange equations are then written as 



- M* - gMrh° ~ 9 P p(r)r 3l ° - e^-^(r) 7 c 



fl> = 0, (10) 



Aa(r) - mla(r) - g 3 a 3 (r) = ^W), (11) 



Aw(r) - mluj(r) - c 3 w 3 (r) = -g u (rh°ii>), (12) 

Ap(r) - m 2 p p(r) = -^(V^r/V), (13) 

AA{r) = _ e $(i^) 7 o^ (14) 

Here, the meson mean fields are functions of r, which is the radial coordinate of the 
nucleon center. The meson mean fields are approximated to be constants over the small 
nucleon volume. We solve the above equations self-consistently with the amount of M* 
and dM*/da obtained in the first step. 



We take the same coupling constants and meson masses as used in nuclear matter to 
study the properties of finite nuclei. We follow some prescriptions in Ref. |25| to deal with 
the center-of-mass corrections and the pair corrections. The charge densities and the rms 
charge radii are calculated by convoluting the point nucleon density with an empirical 
nucleon form factor as in Ref [25|. The calculated results for the binding energies per 
nucleon E/A and the rms charge radii R c are compared with experimental values in Table 
3. For the case of the scalar-vector potential, the calculated results for E/A and R c are 
almost same as the RMF(TMl) results. This is because the derivative dM*/da is almost 
unchanged with the a mean field as shown in Fig. 2. On the other hand, for the case of the 
scalar potential, E/A are somewhat underestimated. This is caused by the increase of g a 
with the a mean field as shown in Fig. 2. We need, instead, the decrease of the effective g a 
in order to achieve good results. In Table 4, the calculated spin-orbit splittings for 40 Ca 
and 20S Pb are presented. The results are similar to those of the RMF(TMl). This is due 
to the fact that the effective mass was used to obtain the coupling constants. In Fig. 6 
we show calculated charge density distributions for Ca in the QMF model and compare 
with the experimental distribution [^7]] and the results in the RMF(TMl) model. Even 
though there are still some discrepancies between the results in the QMF model and the 
experimental values, we consider the QMF model provides reasonable results for finite 
nuclei. 



5 Conclusion 

We have developed the quark mean field (QMF) model to describe the change of nucleon 
properties in nuclei and at the same time the properties of nuclear matter and finite 
nuclei. We have used the constituent quark model for the nucleon, which naturally allows 
the direct coupling of a, u, and p mesons with quarks. The mean field Lagrangian at 
the nucleon level reflects the direct coupling of mesons with quarks merely through the 
appearance of the effective nucleon mass, which is a function of the a mean field. With 
our model setting and our parameter choices, we can perform the numerical calculations 
for nuclear matter and finite nuclei. 

We have investigated the QMF model with different types of the confinement potential. 
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We have taken the mean field Lagrangian with nonlinear terms for both a and uj mesons. 
The comparison between the QMF model and the RMF(TMl) model has been done. The 
QMF model can provides a significantly swollen nucleon radius in nuclear medium. At 
the normal matter density, the nucleon radius increases by about 5% ~ 9%. We calculate 
also the properties of finite nuclei, the resulting binding energies per nucleon E/A and 
charge radii R c are close to the experimental values, as well as the spin-orbit splittings in 
the QMF model are also satisfactory. 
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Figure captions 



Figure 1: The effective nucleon mass M* as functions of the quark mass correction 5m q . 
The results in the QMF model with Xc = \kr 2 are shown by solid curves, while 
those with \c = \kr 2 {\ + 7°)/2 are shown by dashed curves. For each potential 
shown are the two results for two confining strengths. 

Figure 2: The ratios of the u-nucleon coupling in medium, g a ((r) = dM*/da, to that in 
free space, g a (0), as functions of the quark mass correction Sm q . 5m q is connected 
with the a mean field as 5m q = —g%o. The curves are labeled as in Fig.l. 

Figure 3: The energy per nucleon, E/A, as functions of the nuclear matter density p. 
The results in the QMF model with Xc = \kr 2 are shown by solid curves, while 
those with Xc = \kr 2 {\ + 7°)/2 are shown by dashed curves with k = 700. The 
results in the RMF(TMl) model are plotted by dotted curves for comparison. 

Figure 4: The scalar potential, Us, and the vector potential, Uy, as functions of the 
nuclear matter density p. The curves are labeled as in Fig. 3. 

Figure 5: The ratios of the nucleon rms radius R to that in free space Rq as functions 
of the nuclear matter density p. The curves are labeled as in Fig.l. 

Figure 6: The charge density distributions for 40 Ca compared with the experimental 
data (solid curve) ]2^| . The dash-dotted and dashed curves are the results in the 
QMF model with X c = \kr 2 (k = 700) and Xc = \kr 2 (l + 7 °)/2 (k = 700), re- 
spectively. The results in the RMF(TMl) model are plotted by dotted curves for 
comparison. 



14 



Table 1: The nuclear matter properties used to determine the five free parameters in the present 
model. The saturation density and the energy per particle are denoted by p an d E/A, and 
the incompressibility by k, the effective mass by M* and the symmetry energy by a sym . 



Po 


E/A 


k 


M*/M n 




(fm- 3 ) 


(MeV) 


(MeV) 




[MeV) 


0.145 


-16.3 


280 


0.63 


35 



Table 2: The parameters in the QMF model are listed. For comparison, the parameters in the 
RMF(TMl) model are also presented. 



Model 


9% 


9l 


92 


93 


C3 


9 P 


(MeV) 


QMF 


k = 


700 


3.14 


4.20 





50.7 


53.6 


4.3 


470 


Xc = \kr* 


k = 


1000 


2.98 


4.17 





52.8 


36.4 


4.3 


460 


QMF 


k = 


700 


4.18 


4.42 





36.8 


214 


4.3 


515 


Xc = \kr 2 (l + 7 °)/2 


k = 


1000 


4.11 


4.38 





36.4 


167 


4.3 


510 


RMF (TM1) 


g a = 10.029 


9* = 12.614 


-7.2325 


0.6183 


71.308 


4.6322 


511.2 
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Table 3: The binding energies per nucleon E/A and the rms charge radii R c in the present 
model compared with the results in the RMF(TMl) model and the experimental values . 



Model 








E/A 


(MeV) 






R c 


(fm) 










40 Ca 


48 Ca 


90 Zr 


208 Pb 


40 Ca 


48 Ca 


90 Zr 


20S Pb 


QMF 


k = 


700 


7.53 


7.66 


7.92 


7.36 


3.45 


3.46 


4.28 


5.53 


Xc = \kr 2 


k = 


1000 


7.32 


7.46 


7.75 


7.24 


3.45 


3.46 


4.28 


5.53 


QMF 


k = 


700 


8.35 


8.43 


8.54 


7.81 


3.44 


3.46 


4.28 


5.54 


Xc = i£x 2 (l + 7 °)/2 


k = 


1000 


8.21 


8.30 


8.43 


7.73 


3.44 


3.46 


4.27 


5.53 


RMF (TM1) 


8.62 


8.65 


8.71 


7.87 


3.44 


3.45 


4.27 


5.53 


Exp. 


8.55 


8.67 


8.71 


7.87 


3.45 


3.45 


4.26 


5.50 



Table 4: The spin-orbit splittings for Ca and 208 Pb in the present model compared with the 
results in the RMF(TMl) model and the experimental values All quantities are in MeV. 



Model 






A0 Ca 


20S Pb 








Proton 


Neutron 


Proton 


Neutron 








(lcZ 5/2 - W3/2) 


(14/2 - 14/2) 


(I99/2 - I57/2) 


( 2 /7/2 - 2/5/2) 


QMF 


k = 


700 


-5.8 


-5.9 


-3.5 


-1.8 


Xc = \kr 2 


k = 


1000 


-5.8 


-5.8 


-3.5 


-1.8 


QMF 


k = 


700 


-5.6 


-5.6 


-3.3 


-1.9 


Xc = ifcr 2 (l+ 7 °)/2 


k = 


1000 


-5.7 


-5.8 


-3.4 


-1.9 


RMF (TM1) 


-5.7 


-5.7 


-3.4 


-1.8 


Exp. 


-7.2 


-6.3 


-4.0 


-1.8 
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